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Abstract 

Using the renormalon calculus, we study the asymptotic behavior of the pertur- 
bative expansion of the hard-scattering kernels entering the QCD factorization 
formula for the nonleptonic weak decays B'^ D^*^^ AI~ , where M is a light 
meson. In the "large-/?o limit", the kernels are infrared finite and free of end- 
point singularities to all orders of perturbation theory. The leading infrared renor- 
malon singularity corresponding to a power correction of order Aqcd /i^b vanishes 
if the light meson has a symmetric light-cone distribution amplitude. We calcu- 
late the Borel transforms and the corresponding momentum distribution functions 
of the hard-scattering kernels, and resum the series of 0(/3Q~^a") corrections to 
explore the numerical significance of higher-order perturbative and power correc- 
tions. We also derive explicit expressions for the 0(/3o contributions to the 
kernels, and for the renormalon singularities corresponding to power corrections 
of order {Aqcb / n^b)"^ ■ Finally, we study the limit nic — > relevant to charmless 
hadronic decays such as i? — > vrvr. 
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1 Introduction 



The theoretical understanding of nonleptonic weak decays of hadrons is comphcated 
by the intricate effects of strong interactions. Gluon exchange between the quarks is 
characterized by a multitude of relevant mass scales, ranging from the electroweak scale 
/i ~ Mw down to the confinement region /i ~ Aqcd, where perturbative methods fail. 
Recently, however, it was shown that this picture simplifies drastically for most two- 
body hadronic decays of B mesons into final states containing at least one fast, light 
meson 0, 0. In the heavy-quark limit, the decay amplitudes for these processes factorize 
into a semileptonic form factor and a meson decay constant. So-called "nonfactorizable" 
corrections are predominantly perturbative and taken into account by convolutions of 
hard-scattering kernels with light-cone distribution amplitudes of the mesons. Correc- 
tions to this limit are suppressed by powers of Aqcd/'"^6- 

Factorization as established in is a nontrivial property of the decay amplitudes 
that holds true to all orders of perturbation theory and to leading power in AQCD/''^b- 
For practical applications, it is important to obtain an estimate of the leading, power- 
suppressed corrections. Naive dimensional analysis suggests that such corrections are 
of order AQcn/mb ~ 10%; however, it is difficult to quantify this statement. The main 
obstruction is that, in contrast to simpler applications of the heavy-quark expansion to 
exclusive semileptonic decays or inclusive processes 0, the power corrections to factor- 
ization cannot be organized using a local operator product expansion. Whereas certain 
types of potentially large corrections were identified and estimated in |^, it is not (yet) 
possible to write down the complete set of leading power corrections in terms of field- 
theoretic objects. 

The goal of the present work is to gain some insight into the structure of power 
corrections arising from soft, "nonfactorizable" gluon exchange. We use the renormalon 
calculus to investigate the asymptotic behavior of the perturbation series for the hard- 
scattering kernels. These series are divergent and require power corrections of a certain 
pattern in order to be consistently defined. The Borel transforms of the perturbation 
series have singularities in the complex plane, whose structure indicates the power (via 
their position) and strength (via their residues) of these corrections. We study renor- 
malon singularities by adopting the "large-/3o limit" (corresponding to the exchange of 
a single renormalon chain), in which nontrivial partial resummations of perturbation 
series become possible. This method provides the correct location of the renormalon 
singularities, but only approximately accounts for their residues. 

To illustrate our approach, we recall the example of the Adler Z)-function, i.e. the 
Euclidean correlator of two vector currents. For this quantity, the leading infrared (IR) 
renormalon singularity arising from soft gluon exchange indicates a power correction of 
order Aqq^/Q'^, which has the same momentum dependence as the contribution of the 
gluon condensate. A consistent definition of the perturbation series with an accuracy of 
order or better therefore requires inclusion of the gluon condensate, and the value 
of the condensate depends on the resummation prescription (e.g., principal- value Borel 
summation, truncation at the minimal term, etc.) adopted to define the divergent per- 
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turbation series Although the presence of a l/Q'^ power correction proportional to 
the gluon condensate can be inferred from the operator product expansion, it is remark- 
able that the perturbation expansion itself signals the existence of this nonperturbative 
effect through its divergent large-order behavior. The correspondence of renormalon sin- 
gularities and power corrections of (some) higher-dimensional operators in the operator 
product expansion has also been verified with several explicit examples in the context of 
the heavy-quark expansion H, 0, |^ . 

Building on this experience, the renormalon calculus has been applied to observables 
that do not admit an expansion in local operators. Examples are event-shape variables in 

Drell-Yan production [jl3|, |14|, fragmentation func- 



e^e annihilation 



10, 11,0, 13 



tions [T^, [T^ , and structure functions in deep- inelastic scattering [T^, . In these cases 



there is no systematic framework that would allow us to classify the power corrections 
in terms of operator matrix elements. However, the pattern of renormalon singularities 
determines at least a minimal set of power corrections that must be included for a con- 
sistent field-theoretic description. Although, in general, these are not the only sources 
of power-suppressed effects, the inclusion of the corrections corresponding to the lead- 
ing IR renormalons significantly improves the phenomenological predictions. Interesting 
attempts to formalize this method include the dispersive approach developed by Dok- 
shitzer, Marchesini and Webber |[T^ , and the non-local operator method of Korchemsky 
and Sterman . 

Our analysis will be similar in spirit to these approaches in that we will use the 
renormalon calculus to obtain a minimal model of power corrections to factorization in 
hadronic B decays. We will also derive explicit results for the presumably dominant 
part of the two-loop perturbative contributions to the hard-scattering kernels, i.e. the 
contributions of order (3^01- For simplicity, we focus mainly on the class-1 B decays 
into a heavy-light final state such as 5° — * D^*^^M^ ^ where M = vr, p, . . . is a light 
meson. For these processes, the factorization formula takes its simplest form and is best 
established theoretically. 

The effective weak Hamiltonian is (considering h —>■ cud transitions for concreteness) 



K 



cS 



"Ci(/i)Oi(/i) + C2(/i)02(/i) 



(1) 



where 



01 = rfa7^(l - 75)Ma C/37^(l - 75)6/3 ■ 

02 = dajf,{l - 75)^/3 C/37''(l - 75)6a 



(2) 



are local four-quark operators, summed over color indices a, (3. The Wilson coefficients 
Ci{fi) contain short-distance corrections arising from gluon exchange with virtualities 
between the electroweak scale Mw and the renormalization scale fi ~ rrib. The hadronic 
matrix elements of the operators Oi{fi) greatly simplify in the heavy-quark limit mj, ^ 
Aqcd- To leading power in AqcB/^b, and to all orders of perturbation theory, they can 
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be expressed in the factorized form [|I| 
(M-Z}W+|Oi(/i) 15°) 



1 







TiL{x, fi) (M-| J7^(l - 75)m I )(DW+| c7'^(1 - 75)6 



+ T,r{x, fi) (M-| ^7^(1 - j,)u I ) C7^(l + 75)& 



(3) 



$Af(a^5/u) is the leading-twist hght-cone distribution amphtude of the hght meson M, 
normahzed such that /q^ rfx $^(2;, yu) = 1, and x is the longitudinal momentum fraction 
of the d quark in the meson. The current matrix elements on the right-hand side of the 
factorization formula can be expressed in terms of S — D^*^ semileptonic form factors 
and the decay constant /a/ of the meson M. In the heavy-quark limit, "nonfactorizable" 
strong-interaction effects are dominated by hard gluon exchange and included in the 
hard-scattering kernels Ti{x,fi). In these kernels were calculated explicitly at one- 
loop order and were shown to be free of IR divergences at the two-loop order. 

The goal of the present work is to estimate the higher-order corrections in the loop 
expansion of the hard-scattering kernels and, at the same time, to gain some insight into 
the structure of power corrections to factorization. In Section ^ we derive expressions 
for the Borel transforms of the hard-scattering kernels using the single renormalon- chain 
approximation. We present resummation formulae for the kernels in the large-/3o limit 
and establish several important results: the systematic cancellation of IR divergences to 
all orders of perturbation theory, the absence of endpoint singularities in the kernels, and 
the vanishing of the leading renormalon singularity (corresponding to a first-order power 
correction) in the case of a symmetric light-cone distribution amplitude. In Section |^ we 
study higher-order contributions in the perturbative expansion of the kernels. We derive 
exact results for the terms of order Poa^, and give all-order expressions for the anoma- 
lous dimensions of the four-quark operators in the large-/?o limit. We also investigate 
numerically the asymptotic behavior of the perturbation series for the hard-scattering 
kernels. A systematic study of IR renormalons and power corrections is performed in 
Section ^. We present explicit results for the first two IR renormalon singularities cor- 
responding to power corrections of order AqcB/^b and {Aqcb / ■ In Section ^ we 
study the limit rric — > and derive results valid for massless quarks in the final-state. 
They are relevant to charmless decays such as B ^ nn. 



2 Borel transforms and distribution functions 

The perturbation series for the hard-scattering kernels Ti{x,fi) in the factorization for- 
mula d^) can be arranged as 

00 £-1 

T,{x, /i) = + ^ ^ t''\x, /i) /3oX(-") , (4) 

i=l n=0 
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where I is the number of loops, and /?o = y Ca — ^TpUj is the first coefficient of the 
/3-function for an SU(Ai'c) gauge theory {Ca = Nc and Tp = |) with nj hght quark 

flavors. The tree-level coefficients are tf'l = 1, — l/^c, and tf^jl = 0. Our goal is to 
sum the terms of order /3o~^af to all orders of perturbation theory. In other words, we 
consider the limit of large j3o for fixed /3o«s and calculate the kernels Ti{x,fi) to order 
l//3o, neglecting terms of order and higher. Strictly speaking, there is no sensible 
limit of QCD in which (3q may be considered a large parameter (in particular, this is not 
the large- Ai'c limit) - except, perhaps, taking ^/ — > —oo. Nevertheless, retaining only 
the leading terms in l/jSo often gives a good approximation to exact multi-loop results 
(see, e.g., [^), in particular in cases when there is a nearby IR renormalon P^ . 



2.1 Borel summation in the large-/3o limit 

The coefficients of the terms with the highest degree of (3o in (^) are determined by 

diagrams with {£ — 1) light-quark loops, which are rather straightforward to calculate. 
We work in dimensional regularization with 4 — 26 space-time dimensions and adopt 
the MS subtraction scheme. At first order in l/f3o, coupling-constant renormalization is 
accomplished by (/i^ = /i^e'''/47r) 



(47r)' 



6(/i) 



/3oa,(/i2) 



1 



Air 



(5) 



To leading order in the large-/3o limit the "nonfactorizable" vertex corrections to the 
operator Oi in (0) vanish when projected onto color-singlet meson states, i.e. 

T,l{x, /i) = 1 + 0(l//?2) , T,r{x, fi) = 0(l//32) . (6) 

The perturbation series for the kernels corresponding to the operator O2 can be written 



as 



1 



1 + 



2C 



e=i 



\e + b 



(UV subtractions) + 0{1/Pq) 



2Cf ~ FR{e,ie)/ b 



A 



E 



e=i 



e + b 



(7) 



where b = b{fi). The kernel corresponding to the right-handed (i.e. chirality-flipped) 
heavy-quark current in (^) is ultraviolet (UV) finite in the large-/3o limit and does not 
require MS subtractions. 

The functions Fi{e,u) are regular at e = u = 0. Expanding Fi{e,u) in powers of e 
and u, and [b/{e + b)Y in powers of b/e, gives a quadruple sum in (0). Combinatoric 
identities relate the 1/e terms, and hence the anomalous dimensions of the local four- 
quark operators Oi and O2, to the Taylor coefficients of FL{e,0) ||2^. Using the well- 
known fact that the linear combinations 0± = Oi±02 are multiplicatively renormalized. 
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Figure 1: "Nonfactorizable" vertex diagrams contributing to the hard- 
scattering kernels in the factorization formula. The split vertex represents an 
insertion of one of the four-quark operators Oj. The circles represent {i — 1) 
insertions of light-quark loops on the gluon propagators. 



we find that the corresponding anomalous dimensions 7± are given by 

7± = ±^ ^ F,{-b, 0) + 0(1/(3',) . (8) 

The finite terms, which determine the kernels themselves, receive contributions from the 
Taylor coefficients of ^^(e, 0) and Fi{0,u) |2§. We obtain 

2Cp } ^ FL{e,0)-FL{0,0) 



NcT2L{x,fi) = 1 + f de 



-6(m) 





oo 



N,T,nix,i^) = ^ fdue-/'^^^ + 0{l/Pl) . (9) 

un J U 

^ 

If they were unambiguously defined, these expressions would determine the Borel-resum- 
med perturbation series for the hard-scattering kernels in the large-/3o limit. However, 
the presence of singularities located on the integration contour renders the integrals over 
u ambiguous. These (renormalon) singularities will provide us with information about 
power corrections to the factorization formula (Q). 

The functions F^i^e, u) and Fr(£, u) are obtained by computing the "nonfactorizable" 
vertex corrections shown in Figure |I[ It is convenient to introduce the abbreviations 

7T7 — 

z = — , 5={l-z')x, 5={l-z')il-x). (10) 

The variable 5 {5) appears in the calculation of the first (second) diagram in Figure |1[ 
The third and fourth diagrams are obtained using crossing symmetry and require the 
substitutions 

1 5 5 

mb-^rric, z^-, x ^ [1 - x) , 5^ — ^ - , 5-^ — - ie . (11) 

Z Z^ Z'^ 
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These contributions, to which we will refer as "crossed terms" , contain imaginary parts 
(specified by the ie prescriptions), which determine the strong-interaction phases of the 
decay amplitudes. From now on we will omit the ie's; they can be reinstated by recalling 
that z'^ = z'^ — ie. The results are 



FlIe, u) = [ — ] e'' [Die)Y"-^ e, u) - ^2(5, £, u) + crossed terms . 



Fk(£,m)= — e^MD(£)]"/"-^ (5, e, m) + crossed terms, (12) 



where 



Ai{r,e,u) = {1 - ef {r gii + g2o) + u {gio - 920) 

- ur {goo - goi - gio + gn) , 

A2{r,€,u) = (4:- e - e'^){rgn + g2o) +u{gio-2g2o) 

- ur {goo - 901 - 910 + 2gu) , 

Bi{r,e,u) = -ug2o, (13) 



and 



r(m + «-2M)r(i + «) / ^"(1 - v)" ,,,, 

S^n^sUr.e.u)^ r(l+m + n-.-e) /' ^yCr + (1 - ' 

The integral can be expressed in terms of incomplete Euler /3-functions. The notation 
in (|T2|) is such that the functions Ai and Bi correspond to the contributions of the first 
diagram, whereas A2 is obtained from the second diagram. The linear factor of z in front 
of Bi comes from the chirality fiip on the charm-quark line. The contributions from the 
third and fourth diagrams are obtained by performing the substitutions shown in ([TT|). 
Finally, the function 

Di.)^ee^''-^0^^1^ls,Oi.^) (15) 

is related to the contribution of a light-quark loop to the gluon self-energy. 

In the limit u ^ most of the terms in ([13D vanish. However, the functions goo 
and goi have poles at u = 0, corresponding to soft and collinear IR singularities of the 
individual diagrams. Working out the residues we find 



FL{e,0) = e^' [Die)] 



-1 



{3-2e){l + e) ^ 1 x 



crossed terms . (16) 



The logarithm corresponds to a remaining collinear singularity in the sum of the first 
two diagrams, and cancels when adding the crossed terms 0. As a consequence, only 
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the first term in the square brackets, which comes from the first term in the expressions 
for Ai and A2 in (p!3[) , remains. This term depends on the renormalization scheme. 
Specifically, the e-dependent factors in (|13|) follow from the Dirac identities 



l,lula{l - 75) ® - 75) = 4(1 - 7p(l - 75) ® 7^(1 - 75) , 

7^.7.7.(1 - 75) ® 7^V(1 - 75) = 4(4 - £ - e^) 7p(l - 75) ® Yil - 75) , (17) 

which are valid in the so-called "naive dimensional regularization" (NDR) scheme with 
anticommuting 75, and with a projection of evanescent operators as specified in eq. (4.3) 
of The scheme-dependence cancels against that of the Wilson coefficients in the 
effective weak Hamiltonian. Using the expression for D{e) from (p!5|), we obtain 

F(rO)- (3-25)(l + 5)r(4-25) m - 3 hR) 

^^^''°)-"3r(l + 5)P(2-e)r(3-.)' ^M0,0)--3, (18) 



and Ffi{e, 0) = 0. In the opposite limit £ — > we find 



rrib/ 



+ crossed terms . 



Fr(0, u) = f e^"/3 z Bi (S, 0, u) + crossed term . (19) 

We now rewrite the resummed expressions (|^) for the kernels in a more convenient 
form. The /i-dependent factor in the expressions for Fi{0,u) in (|T9|) combines with the 
exponential e'"^^^^^^ in (Of) into the renormalization-scheme invariant combination 



\mbJ \ nib J KnibJ 

where Ay = e^^^Ajj^ is the QCD scale parameter in the so-called "V scheme" [^, and 

ln(m^/A^) 4:71 An 



bv{mb) = , . 2/A2N = JZ = JZ (21) 



is proportional the running coupling constant in that scheme. Here, as always, as(/i^) 
without a label "V" is the coupling constant in the MS scheme. We define the Borel 
transform of the hard-scattering kernel T2/j as 



Sr{u,x) = e 



^ -5u/3 FRjO^u) 



u 



IJ,=mi, U U 



The factor in front of the crossed term comes from the fact that the scale factor 
associated with this term is (/i/mc)^" = -z"^" (/i/mf,)^". With this definition, the Borel 
integral representation for the kernel becomes 



iV,T2^(x,/i) = — ^ /cine-"/''-('"^)5«(n,x) + 0(l//5o'). (23) 



00 



A 



For the kernel T2L, we define the Borel transform as 

-5u/3 

li=mi, 



Sl(u,x) = 



U 



.(12) 



(0,0) 



z 



2u p{34) 



(0,0) 



Ai(5,0,M)-A2(5,0,M) + e-^"/3 1^- 



In- 



1 — X 



+ z 



(24) 

X (crossed terms) 



Here ^^^^^(0,0) 



+ In and f|^^-'(0, 0) = — | — In correspond to the contribu- 



tions of the only first two and the last two diagrams in Figure |l|, respectively, as can be 
seen by taking the limit £ — in (|16|). With the above definition, the two parts of the 
Borel transform corresponding to each pair of diagrams are separately free of UV and 
IR singularities at m = 0, and the result for the kernel takes the form 



N,T2L{x,^l) = 1 + 



2Ch 



2Cf 
2Cf 



3 In 



de 

e 

b{mb) 



(3-2£)(l + e)r(4-2e) 

3r(i + £) r2(2-£)r(3-e) 



3 , X 

- + In 

2 1-x 



In 



b{mb) 



du e-"/''^(™'') Sl{u, x) + 0{l/(3l) . 



(25) 



In the above expressions for the hard-scattering kernels, the Borel integrals (the 
integrals over u) contain all nontrivial information about the asymptotic divergence of 
the perturbation expansions. They are defined in a renormalization-scheme invariant 
way. So-called IR renormalon singularities located on the integration contour along the 
positive real w-axis render these integrals ill-defined, and hence the perturbation series 

In the large-/3o limit, renormalon singularities show 



are not Borel summable ||2^, ^ ^ | . 
up as poles at half-integer values of u, as seen from the explicit form of the functions 
gmn in (0)- A pole singularity located at u = k/2 corresponds to a factorial growth 



tf ~ {k/2)-^ (^ - 1)! of the perturbative expansion coefficients in (|^), and leads to 
an irreducible ambiguity of order [Aqcd /tTT'b)^ in the definition of the resummed series. 
From the expression for Qmn it follows that for n = these functions have single poles 
at M = i. Working out the residues we obtain 



Sl{u,x) 



4 



1 



1 -2m 



+ regular terms . 



(26) 



whereas Sr{u, x) is regular at n = i. A pole at u = ^ corresponds to a power correction of 
first order in AqcB/^b, which could be large for realistic heavy-quark masses. However, 
because the residue of the pole in ( pUf ) is antisymmetric under exchange of x ^ (1 — x), 
this contribution vanishes for the case of a symmetric light-cone distribution amplitude. 
Therefore, the renormalon analysis does not indicate a first-order power correction for 
decays such as — *• D^*^^n^ or i?*^ D^*^^ p~ . A more detailed investigation of power 
corrections will be presented in Section ^. 
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2.2 Momentum distribution functions 



The analysis of higher-order renormalon singularities and the numerical evaluation of 
the Borel integrals in the form of ( p3l) and (^) are difficult because of the complicated 
structure of the Borel transforms Si{u,x). It is advantageous to rewrite the Borel inte- 
grals as integrals over a running coupling constant multiplied by weight functions, using 
the formalism developed in [p^] : 



dr 



aATe 



-5/3 



mt 



4:71 



(27) 



This elucidates the connection between IR (UV) renormalons and small-momentum 
(large-momentum) contributions in Feynman diagrams. The functions Wj(r, x) deter- 
mine the distribution of the gluon virtualities inside the vertex-correction diagrams in 
Figure |l|. Technically, they are the inverse Mellin transforms of the Borel images Si{u, x). 
The running coupling constant under the integral on the right-hand side is nothing but 
the coupling a^^^rml) in the V scheme, as defined in (|21|). Equation (pTf ) has been 
established in the literature as an integral representation of the Borel sum for quantities 
defined in Euclidean kinematics It is applied here for the first time to a situation 



in which the quantity of interest has a more complicated analytic structure and a non- 
vanishing dispersive part. We have checked numerically that the two sides in (^71) agree 
even if the imaginary part is nonzero. 

The distribution functions Wi{T, x) can be computed in terms of Feynman parameter 
integrals starting from the relation [p2| 



/(It 
— Wi{T,x) r"", 





which is valid for — 1 < Rew < The result is 



(28) 



WL{r, x) = fUr, 5, 5) + hir/z', -5/z\ -5/z') 



wr{t, x) = z fn^r, 6) + - fR^T/ z^, 



(29) 



where the first (second) terms on the right-hand side correspond to the first (last) two 
diagrams in Figure ^ We find 



+ 



25 5y 5 
V 



6^ 



6 



- 1 



{6^6} 



fRir,S) = T 



(30) 
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where 9{x) is the step function, and 



Note that the terms inside the square brackets in the expression for fi are antisymmetric 
in X (1 — x) and thus vanish for the case of a symmetric hght-cone distribution 
amphtude. 

The above expressions for the distribution functions Wi{T,x) are a central result of 
this work. They will allow us to extract information about renormalon singularities (and 
hence the structure of power corrections) as well as about terms of arbitrarily high order 
in perturbation theory. Let us now stress some important general properties of these 
functions. 

IR cancellations: Expanding the distribution functions for small r shows that they 
vanish at least as fast as ^/T as r — 0, and therefore the integrals in (|27| ) are convergent 
in the IR region. As a consequence, the resummed expressions for the hard-scattering 
kernels in ( p3D and (|25|) are IR finite. This is a nontrivial result of our analysis, which 
demonstrates the IR finiteness of the hard-scattering kernels (and hence factorization) 
to all orders of perturbation theory in the large-/5o limit. (After MS subtractions, ac- 
complished by the term proportional to 6'(re^/^ — 1) in the expression for wl{t,x), the 
distribution functions vanish like 1/r for large r, and hence the integrals are also UV 
convergent.) 

Endpoint behavior: The proof of IR finiteness of the hard-scattering kernels alone 
does not establish factorization; in addition, one must show that the convolutions of 
the hard-scattering kernels with the light-cone distribution amplitude of the meson M 
are convergent. In |^, an explicit calculation showed that at one-loop order the hard- 
scattering kernels tend to a constant (modulo integrable endpoint logarithms) for x — >■ 
or X ^ 1. We find that this property persists for the resummed perturbation series in 
the large-/3o limit. Therefore, the integrals over the light-cone distribution amplitude 
^m{x) are convergent at the endpoints. (Endpoint singularities at least as strong as 
l/x^ or 1/(1 — would be required to spoil factorization.) We define 

1 

Wl{t) = Jdx $A/(x) wl{t, x) = Fl{t, 1 - ^2) + Fl{t/z'', 1 - 1/^2) , 



} 1 

WRiT)= dx^M{x)wR{r,x) = zFr{t,1-z^) + -Fr{t/z\1-1/z^), (32) 
J z 



where the functions Fi are the convolutions of fi with the light-cone distribution am- 
plitude. As an example, we perform these convolutions adopting the asymptotic form 
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Figure 2: Integrated distribution functions Wifi^r) obtained using the asymp- 
totic light-cone distribution amphtude. Sohd (dashed) hues show the real 
(imaginary) parts. 



^o{x) = 6x(l — x) of the light-cone distribution amplitude, valid for light pseudoscalar 
and vector mesons in the limit ^ oo. This yields 

^ / ,x 3r/„ 12r7\ 6r / r(rf + 2) 3/7 \ , / r]d\ 

+ 5l!(l±^L,(-!^)+|.(re-/»-l). 
. |(_, , „ _ , , ^) , (, , i^) _ 5^u.(- i^?) , (33, 

where Li2(x) = — jQ{dt/t) ln(l —t) is the dilogarithm, and the subscript "0" refers to the 
asymptotic distribution amplitude. In Figure @, we show the real and imaginary parts 
of the functions Wi^lr) for the case where z = rric/mb = 0.3. Note that the steps in the 
real part of the left-handed kernel are an artifact of the MS subtractions applied to the 
Borel transforms. 



Strong-interaction phases: The hard-scattering kernels contain imaginary parts due 
to gluon exchange among the final-state quarks in the third and fourth diagrams in 
Figure 0. These imaginary parts determine the strong-interaction phases of the decay 
amplitudes in the heavy-quark limit. They are obtained from the branch cuts of the 
logarithms in the crossed terms in (|29|). Imaginary parts exist for r < 5^/(z^ + 6) or 
r < 5^/ {z"^ + S), which in any case implies that r < (1 — z^y. 

Perturbative expansion: If the running coupling constant under the integral in (^) 
is expanded in powers of a fixed coupling constant using 

a.ire-^/^ml) = a.(/.^) £ - log ^ - lnr)"(^^^)" , (34) 
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the integral reproduces the perturbative expansion in powers of Q;s(/i^). Computing 
the first n terms of the perturbation series in the large-/3o hmit requires evaluating the 
integrals /g°° (dr/r) (In r)"~^Wi (r, x) over the distribution functions. To perform these 
integrals analytically, it is convenient to change variables from r to rj. Below, we derive 
explicit expressions for the hard-scattering kernels at order /9o«s- 

The shapes of the distribution functions determine the momentum regions from which 
there arise important contributions to the Borel integrals. As a first indicator, one 



may consider the average value of In r, which determines the so-called ELM scale p6 
through /uIlm = e^'"^"^^ m^. (This relation holds in the V scheme. In the MS scheme, 
the BLM scale /xblm is reduced by a factor e~^/^.) As explained in p^, a BLM scale is 



meaningful only in cases where the quantity of interest is renormalization-group invariant 
and has a distribution function of definite sign. In our case, this means that we can 
introduce BLM scales only for the imaginary parts of the kernels. The dashed curves 
in Figure ^ show that these scales are significantly smaller than mi,. Specifically, we 
find /iBLM ~ 0.33mf, (0.39m;,) for the imaginary part of the left-handed (right-handed) 
kernel. When transformed into the MS scheme, the corresponding scales are of order 
0.6-0.7 GeV. This suggests that the imaginary parts of the kernels, and hence the strong- 
interaction phases of the decay amplitudes, are susceptible to nonperturbative physics. 

Renormalon ambiguities and power corrections: The integrals over the distribu- 
tion functions in (p7|) are ambiguous because of the Landau pole in the running coupling 
constant, 

47r lnr + ln(m2/A2.) ' ^'^^> 

which is located at Inr = Inr^ = ln(AY/m^). The principal values of these integrals 



exactly reproduce the principal values of the original Borel integrals The residue 
of the Landau pole therefore provides a measure of the renormalon ambiguity, which we 
define as 



2Cf . . A2 
— — u;,(ri,x) ; = - 
Po 



\NcT2i{x)\ = ——Wi{TL,x) ; tl = — . (36) 



Since <^ 1, the leading contributions to the ambiguity are given by the first few terms 
in the Taylor expansion of the distribution functions for small r. 



3 Large-order behavior of perturbation theory 

In this section we study in more detail the significance of higher-order perturbative 
corrections to the hard-scattering kernels, using the all-order results derived in the large- 
limit. We start by focusing on the anomalous dimensions of the four-quark operators 
Oi and O2, whose perturbation series are convergent and can be summed exactly in 
the large-/3o limit. We then consider higher-order contributions to the hard-scattering 
kernels. 
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3.1 All-order results for the anomalous dimensions 



Using the exphcit result for the function File, 0) in (0), we obtain from 

^NcTl b 2(3 + 26)(l-6)r(4 + 26) 



7± 



N, /?o 3r(l-6)r2(2 + 6)r(3 + 6) 



+ 0(1/(3^) 



6 



Air 



Atx J 



6 



47r/ 



+ . . . 



(37) 



The radius of convergence of the perturbation series is /?o|cts| < 47r. The all-order results 
t may 1 



in the large-/3o limit may be compared with the exact two-loop expressions 

2 



7± 



a. 



« - , ^ ^ , , 57 21\ 



An) 



+ 



(38) 



Numerically, keeping only the term proportional to (3q a1 is an excellent approximation 
for 7_|_ but not for 7_. However, in both cases the two- loop coefficients are small, so there 
is no reason to expect a dominance of the /3o terms. (These terms are often dominant in 
cases where the series is divergent and the expansion coefficients are large.) 



3.2 Partial two-loop results for the hard-scattering kernels 



Perturbative expansions of the hard-scattering kernels can be obtained by expanding the 
running coupling constant under the integral in (P7|), as well as the couplings asijnl) 
and as{m?^ in (p5|), in powers of ^^(/i^). The kernels are then written as 



T2.(x,/.) = i-5.. + ^ 



An 



t«(x,/.)+(^)>otf(x, /.) + ... 



(39) 



where i = L,R, and 6iL = 1 if i = L and zero otherwise. The results can be expressed 
in a compact form by introducing the functions 



1 



+ 



In 5 

1-5 ' 

5\xi5 



\n^6 - Li2(l -5) -{5 ^5] 



2(1-5) 2(1 -5)^ 



5 



71n(5-61n25-6Li2(l -5) 



4(1-5) 

\W6-\ti5 , ,,,2. 2,3^ Li2(l-(5) 

+ — r- + ln(l - 5) \tl^5 - - \t?5 + — 

1—0 3 1—0 



+ Li3(l-(5) + 2Li3((5)-{5^5} 
3 5 



+ 



4(1-5) 4(1-5)2 
13 



ln5-21n^5-2Li2(l -5) 



(40) 



Here Li3(x) = jQ{dt/t)Li2{t) is the trilogarithm function. At one-loop order we repro- 
duce the expressions obtained in [||], i.e. 



t'i\x,i2) = -9 + 61n — + (3 + 21n 



nib 



X 



1 — X 



In z + h1\6, 5) + h^l'{-5/z\-5/z 



(1), 



t^^{x,^^)=zh'^>{5) + -h'^\-5/z'). 



(41) 



For the corrections of order f3o we find 



t^l\x,^i) 



395 



24 



-7r^-61n^ 



+ 171n 



3 + 2 In ■ 



X 



1 — X 



In 2; + 2 In 2; In 



nib 



+ 



17 10 , X 

\ In 

2 3 1-x 



\nz 



--2\n^)hi\S,5)+(--2\n^-2\nz) h^^ {-5/z\ -5 / z 



hf{5,5) + hf{-5/z\-5/z') 



,(2) 



- - 2 In — 

3 



zh'~^{b)^-h'^{-blz') 



nib 



/i 



I - 2 In — - 2 In ;z ) ^ /i^J V V^^) 



,(2), 



(42) 



395 

■ 24 



in the expressions for t^l\x,^) and t^^'{x,jj) are scheme 



(2), 



The constants —9 and 
dependent and specific to the NDR scheme. 

In Figure ^ we show the one-and two-loop kernels for z = 0.3 and two values of 
the renormalization scale. In order to suppress the integrable, logarithmic endpoint 
singularities we have multiplied the kernels with the asymptotic distribution amplitude 
'^'0(3;) = 6x(l — x). The two-loop contributions are sizable. They are suppressed rela- 
tive to the one-loop contributions by a factor /?o«s/4vr ~ 0.15 (for fj, ^ nib), but have 
coefficients that are typically larger by about a factor 3. 



3.3 Numerical results 

The results derived in Section ^ allow us to study the significance of higher-order pertur- 
bative contributions to the hard-scattering kernels in the large-/3o limit, which may be 
considered as a toy model for the asymptotic behavior of perturbation theory. In Table |I], 
we show the results for the kernels obtained at fixed order in perturbation theory and 
compare them with the principal values of the Borel-resummed series computed from 
(p3|) , (p5|) and (^) . For simplicity, the kernels are integrated over x with the asymptotic 
light-cone distribution amplitude. As a measure of the irreducible uncertainty in the 
value of the Borel integral we quote the renormalon ambiguity defined in (p6|). As input 
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Figure 3: One-loop (upper plots) and partial two-loop (lower plots) contribu- 
tions to the hard-scattering kernels for /j, = (black curves) and, if different, 
H = mh/2 (gray curves). Results are multiplied by the asymptotic light-cone 
distribution amplitude. Solid (dashed) lines show the real (imaginary) parts. 

parameters we take nib = 4.2 GeV, z = rric/mb = 0.3, and Uf = 4. We use the one-loop 
running coupling constant (as is appropriate in the large-/3o limit) with = 150 MeV, 
which gives ^^(m^) = 0.226. The corresponding scale parameter in the V scheme is 
Ay = 345 MeV. 

The results in the table show that the large-order behavior of the series is improved 
by lowering the renormalization scale. For /x = nib/ 2, the two- and three-loop results 
(A^ = 2 or 3) are reasonably close to the Borel-resummed value of the entire series. (We 
stress that these resummed values must be considered as estimates only. As a rule of 
thumb, the "true" resummed values lie somewhere between the two-loop results and the 
resummed values obtained in the large- /Sq limit.) An exception is the real part of the 
right-handed kernel, which suffers from large cancellations in the integral over the weight 
function. In this case, higher-order contributions in both ag and l/(3o are potentially 
more important. Note that the imaginary parts of the kernels approach their asymptotic 
values slower than the real parts. This is a reflection of the low BLM scales obtained 
earlier. Physically, it means that the strong-interaction phases of the decay amplitudes 
are sensitive to nonperturbative hadronization effects. 
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Table 1: Fixed-order 0{a^) perturbative approximations to the hard-scattering kernels 
(convoluted with the asymptotic distribution amplitude) and the corresponding principal 
values of the Borel-resummed series, for two choices of the renormalization scale. Results 
for the right-handed kernel are multiplied by 10. 





J^dxMx)N,T2L 




lOx J^dx%{x)N,T2R{x,ii) 


N 


11^ nib 




mb/2 




II = ■mb/2 





1 




1 








1 


0.53 - 0.18i 


0.66 


- 0.23i 


-0.05 - 0.37i 


-0.06 - 0A7i 


2 


0.39 - 0.29i 


0.54 


- OMi 


0.02 - 0.57i 


0.06 - 0.66i 


3 


0.34 - 0.35i 


0.52 


- OAOi 


0.09 - 0.68i 


0.15 - 0.76i 


4 


0.32 - 0.40i 


0.52 


- 0.45i 


0.17 - 0.75i 


0.25 - O.Sli 





0.31 - O.Ui 


0.52 


— OAm 


0.24 - 0.80i 


0.33 - 0.85i 


6 


0.31 - 0.48i 


0.52 


- 0.53i 


0.31 - 0.83i 


0.43 - 0.88i 


7 


0.31 - O.bli 


0.53 


- 0.58i 


0.39 - OMi 


0.55 - 0.90i 


8 


0.32 - 0.55i 


0.54 


- 0.65i 


0.49 - 0.88i 


0.73 - 0.93i 


9 


0.33 - 0.60i 


0.58 


- 0.76i 


0.61 - 0.91i 


0.99 - 0.96i 


10 


0.35 - 0.67i 


0.62 


- OMi 


0.78 - 0.93i 


1.46 - OMi 


Borel Sum 


0.29 - 0.47i 


0.49 


- 0.47i 


0.26 - 0.90i 


0.26 - 0.90i 


Aiiil)iguit}' 


0.003 - 0.01/ 


0.003 


- 0.01/ 


0.10 - 0.03/ 


0.10 - 0.03?: 



Another important observation from the table is that the renormalon ambiguity, 
given in the last row, is always much smaller than the contributions of the first few 
terms in the perturbation series. (Note that the very small ambiguity for the real part 
of the left-handed kernel results from cancellations between various terms in the Taylor 
expansion of the function Wl(tz:). The leading power contribution alone is 0.017.) The 
series for the kernels are well-behaved in the sense that their divergent behavior sets 
in only in high orders {N ~ 6). This suggests that power corrections to factorization 
(at least those connected with "nonfactorizablc" soft gluon exchange) arc smaller than 
perturbative contributions due to hard gluon exchange. For this reason, it would be 
useful to have complete two-loop expressions for the hard-scattering kernels. 

As a final remark, we note that for the real part of the left-handed kernel the perturba- 
tion series never comes close (to within an amount of order the renormalon ambiguity) to 
the principal value of the Borel integral. This is because the kernels relevant to 6 ^ cud 
transitions are the sum of two perturbation series, one with a "natural scale" of order 
nib (diagrams 1 and 2), and the other with a "natural scale" of order nic (diagrams 3 
and 4). These two series reach their minimal term at different values of A^. When their 
sum is truncated at fixed N, one of the two series can be close to its minimal term, but 
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not both. This is a general feature of multi-scale problems in quantum field theory. If 
the scales ruf, and were widely separated, one could disentangle the two series by 
constructing an effective field theory. However, this does not help if, as in the real world, 
the two scales are relatively close to each other. 



4 IR renormalons and power corrections 

After the study of higher-order perturbative contributions in the previous section, we now 
turn to the discussion of power-suppressed effects. We stress, from the outset, that the 
power corrections inferred from the renormalon analysis arise from soft "nonfactorizable" 
gluon exchange of the type shown in Figure ||. These are not the only sources of power 
corrections to the factorization formula (|^). For instance, different decay topologies such 
as weak annihilation or gluon exchange with the spectator quark in the B meson are 
known to contribute at order AqcD/w^fe and have been estimated in For the case 
of D^*^~^M~ decays, their effects were found to be less than 10%, as suggested 

by naive power counting. (For the decay B —>■ tttt, the leading power corrections to the 
factorization formula have recently been estimated using light-cone QCD sum rules. The 
result is, again, a moderate correction of less than 10% PT| .) 

As mentioned earlier, due to the presence of IR renormalon singularities the re- 
summed perturbation series for the hard-scattering kernels can only be defined up to 
an irreducible ambiguity, which in the large-/3o limit can be estimated from (|36|). We 
will now analyze the structure of the leading power-suppressed effects in more detail, 
and obtain a minimal model of power corrections that is consistent with the renormalon 
analysis. It follows from (|36|) that the renormalon ambiguity is given in terms of the dis- 
tribution functions Wi{T, x) evaluated at the small value t = = (Av/mt)^ ^ 1. Terms 
of order r'^/^ in the Taylor expansions of these functions correspond to power corrections 
of order (Aqcd/'"^6)'^- We focus first on the functions /j in (^) and note that the limits 
r ^ and x — > 0, 1 do not commute. Naive expansions of these functions for r 
would yield more and more powers of 1/x and 1/(1 — x), thereby introducing endpoint 
singularities. However, we have seen earlier that the kernels can be integrated over x 
without encountering such singularities. Therefore, the expansions for small r must be 
written in terms of distributions. For simplicity, we shall assume that the light-cone 
distribution amplitude ^m{x) vanishes linearly at the endpoints. This is true for the 
asymptotic distribution amplitude and, more generally, whenever the amplitude can be 
approximated by a finite expansion in Gegenbauer polynomials, i.e. 



$M(a;,yu) = Qx{l - x) 



N 



1 + E «nV) CPHSX - 1) 



n>l 



(43) 



The Gegenbauer moments a^{n) are multiplicatively renormalized and vanish for /i 
CO, so that $Af(a^, At) tends to its asymptotic form. Here we adopt the Gegenbauer 
expansion for convenience only; our analysis of the endpoint behavior in Section ^ showed 
that the kernels can be integrated with any normalizable and smooth function ^m{x)- 
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In our analysis we keep only the leading terms of order ^JT or r, neglecting higher- 
order contributions corresponding to power corrections of order (Aqco/'^^fe)^ and higher. 
We obtain 



ax\ ax 2 



' _^A(x,v^/ci)-{x^(l-a;)} 



dx d?x 



fR{r,6) 



T 



2dx 



+ 0{r 



3/2 N 



where d 



z^, and we have introduced the distribution 



A(x, a) 



X 



— In a; — In a + - 



+ 



5{x) 



, 2 , 1 
in a-lna + — + - 



Its integral with the light-cone distribution amplitude is defined as 



(44) 



(45) 



j dx ^^^^^ A(x, a) = j dx 



^m{x) 



X 



- *m(0) 



In X — In a + 



+ 



*m(0) 



In a — In a 



vr^ 1 
Y+2 



where 



*m(0) = lim 



a:->0 X 



6 



n>l 



(46) 



(47) 



The -distribution for x ^ 1 is defined in a similar way. Its evaluation involves the 
derivative $^(1), which up to an overall sign is given by the same sum over Gegenbauer 
moments, but without the factor (—1)"". 

To get the corresponding expansions of the functions Wj(r, x) we must add the crossed 
terms, for which t ^ t/z^ and d — > —d/z^. This yields 



wl{t,x) 



2r 



dil — X 



(\\iz + in) 



+ 



2v^ 



'1 — z) X d^x 



A(x, ^/r/d) - A(x, -zy/r/d) 



dx 



x)} 



+ 0(r3/2) 



Wr{t,x) 



T / Z 



2d \x z{l — x) 



(48) 



The asymmetric part of the function wl has a small-r behavior corresponding to a first- 
order power correction in Aqcb /rnb, in accordance with the result (|2^) for the leading 
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renormalon pole. When integrated with a symmetric hght-cone distribution amphtude 
this term vanishes, and the leading power corrections are of order (AgcD/'^^fe)^- 

In the limit of isospin symmetry, the light-cone distribution amplitude for a pion or 
a p meson is symmetric in x ^ (1 — x). We then obtain 

W^L(r^) = ^^fln^-m)fl+ af)+0[(AQCD/m,f], 



m;; — ml \ mr. 



Wn{r,) = ( 1 + E + 0[(AqcdM)3] , (49) 



2mbmc 



n=even 



where we have used that jQ{dx/x) $M(a^) = 3(1 + + a^^ + . . .) for a symmetric 
distribution amplitude. For decays involving strange particles in the final state, such as 
Bq — s> K~ , the function ^m{x) is no longer expected to be symmetric. In such a 
case, the function Wl{t) receives a first-order power correction given by 

12Av / .r-^ M 



W^tl) = E ^n] + 0[(AqcdM)2] . (50) 

A comment is in order concerning the peculiar dependence of the power corrections 
in ( ^9D and (0) on the heavy-quark masses, which prevents us from taking the limits 
mc — > mb or m-c — > 0. The limit where the charm-quark mass tends to zero is actually 
not a singular one, but has to be dealt with carefully. We will come back to this in 
the following section. The fact that some power corrections blow up for m^ — *■ m^ is 
physical. There are three relevant mass scales in this problem: the mass of the decaying 
h quark, the mass of the charm quark, and the energy Em ~ {fnl — m^)/(2mb) of the 
light final-state meson. The factorization properties of the decay amplitude in the heavy- 
quark limit rely crucially on the fact that Em S> Aqcd, since only then can the color 



transparency argument [^2|, ^ be employed to demonstrate the cancellation of soft IR 
contributions [Q]. This condition no longer holds in the limit m^. m^. 

Finally, we can use the above results to write down a minimal model for the power 
corrections due to soft gluon exchange in hadronic B decays. The decay amplitudes for 
the class- 1 decays B'^ — > D^*^^M~ are conveniently parameterized in terms of quantities 
ai(D*^*)M), which contain the QCD corrections to the results obtained using "naive" 
factorization. The factorization formula (^) implies that 

1 

ai = Y Ci{p) dx^M{x,iJ,) TiL{x,p)±TiR{x,p) + ©(AqcoMfe) , (51) 

where the upper (lower) sign refers to the case with a D {D*) meson in the final state. 
Inserting the Borel-resummed expressions for the hard-scattering kernels, we find that 
the renormalon ambiguity in ai is given by 



Wl{tl)±Wr{tl) +0(1//32). (52) 
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To obtain a model for the leading power corrections, we insert here expressions (^) 
or (pO|), depending on whether or not the distribution amplitude is symmetric, and 
replace Ay by nonperturbative hadronic parameters Al,r- Ideally, these parameters 
would be determined from a fit to experimental data. However, for A^^/j ~ 0.5 GeV and 
a symmetric light-cone distribution amplitude, we find power corrections of order few 
times 10^^, which are insignificantly small. This is in part due to the smallness of the 
Wilson coefficient C2. The power corrections could be much larger for other decays, such 
as i? ^ TTTT. For decays with a strange meson in the final state, such as — > D~^K~, 
our model predicts first-order power corrections as large as a few percent. Unfortunately, 
these modes are Cabibbo suppressed, and no experimental data has been reported yet. 



5 Results for charmless hadronic decays 

Our discussion so far has referred to the simplest application of the QCD factorization 
formula, i.e. to class- 1 B decays into a heavy-light final state. However, factorization 
in the heavy-quark limit also occurs for the phenomenologically more interesting, rare 
hadronic decays into two light mesons. Examples are the charmless decays S — > tttt 
and B —>■ ttK, which might provide information about the CP-violating phase of the 
quark mixing matrix [§. The effective weak Hamiltonian for these processes contains 
many penguin operators besides the current-current operators Oi and O2, and several 
different decay topologies must be considered in addition to the diagrams shown in 
Figure |I|. A complete renormalon analysis for these processes is beyond the scope of 
this paper; however, the "nonfactorizable" vertex corrections investigated here are an 
important part of such an analysis. It is thus interesting to apply our results to the case 
where the charm quark in the final state is replaced by a massless u quark. 

The limit z — >■ is smooth but must be taken carefully. For instance, the power 
corrections in (^) cannot be extrapolated to rric — >■ 0, but we obtain regular expressions 
by first computing the distribution functions in the limit z ^ 0, and then expanding the 
results for small r. We now collect the most important formulae valid for z = 0. 



5.1 Borel representation and distribution function 

In the massless limit the chirality of the external quark states is preserved, and hence 
the right-handed kernel Tr vanishes for z 0. This can also be seen explicitly by taking 
the limit z — > in our results for this kernel. The Borel-resummed expression (|25|) for 
the left-handed kernel simplifies to 



N,T2Lix,fJ.) = 1 + 



2Cf 



de 

e 



(3-25)(l + £)r(4-2£) 
3T{l + e) T^{2-e)T{3-e) 



+ 3 1n-^^ 

b{mb) J 



„ r dr , 
2Cf J —wl{t,x^ 



47r 



(53) 
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Figure 4: Integrated distribution function Wl{t) for rric = 0. The left-hand 
plot shows the asymptotic result, and the right-hand one the contribution 
proportional to the first Gegenbauer moment. 



where 



WL{r,x) 



2x 



x-^ 



+ 



x^ 

T — 1] 



X 



le] +3^(re-^/3 - 1) 



1-^ In 1 + ^ + l-I In 

X \ X'^ J \ T J \ Xj \ T 



— {x ^ (1 — x)} 



(54) 



As before, the kernel is well-behaved in the endpoint regions a; ^ and x —>■ 1. In 
particular, the integral of wl{t, x) with the asymptotic light-cone distribution amplitude 
yields 

WL,o{r) = (1 - 13r]) - 6r(l + 3r] + 3r) In (l + ^) + 6r(l - r) In (^1 - ^ - ze 



+ 12t' Li2 (- ^) + Li2 + ^e) + 3 0(r e''^' - 1) 



(55) 



This function is shown in the left-hand plot in Figure ^. In the second plot, we show 
the contribution proportional to the first Gegenbauer moment, denoted by Wl,i{t). It 
is evident that the real part of Wl,i{t) receives large contributions from the region of 
very small r, corresponding to low gluon virtualities. This is in accordance with the fact 
that the renormalon ambiguity for this contribution is of first order in Aqco/'^fc- 
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5.2 Partial two-loop results for the hard-scattering kernel 



The easiest way to obtain the one and two-loop coefficients of the kernel in the massless 
case is to take the limit 2; — > in the expressions collected in Section |3.2| . It is a nontrivial 
check of our results that all logarithms of the mass ratio z cancel in that limit. For the 
expansion coefficients in ( ^9|) we obtain 



t^]j) {x, /u) = 6 In 



rrib 



-9 + 61n— + ^ 
fi 2 



3/1-22; 



1 — X 



In X — iiT 



+ 



In^x 



2 



+ 



In a; 

1 — X 



+ Li2(x) — [ - + in \ \nx — {x ^ {1 — x)} 



llln^-^+f^ 
24 V3 



2ln^)ti\x,f.) 



3(1 -3a;) , 2 7(1 -2a;) , 3(1 - x) 



+ 



4(1 -x) 
" In^a; 

l-3a; 



4(1 -x) 



2x 



Li2(x) +i7r[- Inx - - 



In^x ln(l — x) 



1 + 3x 
4(1 -x) 



In X + TT 



Inx 



ITT 



Li2(x) - Li3(x) + — (In x - 31nx) - {x ^ 1 



x} 



(56) 



5.3 Numerical results 

In Table we show the results for the kernel obtained at fixed order in perturbation 
theory and compare them with the principal value of the Borel-resummed series. To 
illustrate the effect of the leading IR renormalon singularity, we include the contribution 
of the ffist Gegenbauer moment a'^ . Then the light-cone distribution amplitude is no 
longer symmetric, and the leading renormalon ambiguity is of ffist order in AqcD/'^ft- 

The results obtained using the asymptotic distribution amplitude exhibit a similar 
behavior as in the case with a finite charm-quark mass. The convergence of the series 
is improved by using the smaller renormalization scale /i = mf,/2. In this case, the two- 
loop result is already very close to the asymptotic value. As previously, the imaginary 
part of the kernel reaches its asymptotic value at larger than the real part. The 
situation for the contribution proportional to the ffist Gegenbauer moment is different. 
Whereas the imaginary part shows a similar behavior as for the leading term, the series 
for the real part of Wl,i{t) diverges already in low orders. Even the two-loop results 
exceed the asymptotic values, and starting at A^ ~ 4 the expansion coefficients exhibit 
a rapid factorial growth. This is a reflection of the leading renormalon pole at m = |, 
corresponding to a ffist-order power correction to the real part, which is absent in the 
case of a symmetric distribution amplitude. Accordingly, the renormalon ambiguity is 
an order of magnitude larger than in the symmetric case. 

Finally, we note that using the one-loop expressions for the kernel, as is done in 
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Table 2: Fixed-order 0{a^) perturbative approximations to the hard-scattering kernel 
in the limit ttIc —>■ and the corresponding principal value of the Borel-resummed se- 
ries, for two choices of the renormalization scale. Shown are the first two terms in the 
Gegenbauer expansion of the distribution amplitude. 











X,fl) 








fi = nih 


fi = mb/2 









1 






1 






1 


(0.56 


-0.23i; 


+ af (0.13 


-0.23i) 


(0.69- 


0.292) + af 


;0.17- 


0.29i) 


2 


(0.44 


-0.35f 


+ af (0.31 


-0.35i) 


(0.61- 


0.41z) + af 


;0.41 - 


OAOi) 


3 


(0.41 


- 0A3i] 


+ af (0.51 


- 0.42z) 


(0.61 - 


0.48z) + af 


^0.67- 


0.46z) 


4 


(0.41 


- 0A8i] 


+ af (0.77 


- o.m) 


(0.62- 


0.53i) + af 


;i.06- 


OMi) 





(0.42 


-0.52i; 


+ af (1.13 


-0.50i) 


(0.65- 


0.57i) + af 


;i.71 - 


OMi) 


6 


(0.44 


-0.55f 


+ af (1.76 


-0.53i) 


(0.67- 


OMi) + af 


^3.05- 


OMi) 


7 


(0.46 


-0.59i; 


+ af (2.96 


-0.57i) 


(0.72- 


OMi) + af 


^6.28- 


OMi) 


8 


(0.49 


-OMi] 


+ af (5.63 


-0.60z) 


(0.78- 


0.72z) + af 


;i5.2 - 


OMi) 


9 


(0.54 


-OMi] 


+ af (12.3 


-0.65i) 


(0.91- 


0.82z) + af 


^42.7- 


0.79z) 


10 


(0.61 


- 0.74f 


+ af (30.5 


-0.71i) 


(1.15- 


0.99i) + af 


;i38. - 


0.96i) 


Borel Sum 


(0.38 


-0.55f 


+ af (0.27 


-0.53i) 


(0.59- 


0.55z) + af 


;0.27- 


0.53z) 


Ambiguity 


(0.02 


- 0.04i; 


+ af (0.19 


- OMi) 


(0.02- 


OMi) + af 


;0.19- 


OMi) 



all phenomenological applications of the QCD factorization approach to date, gives a 
reasonable approximation to the real part (for /i = mb/2), but underestimates the strong- 
interaction phase by almost a factor 2. This observation may be relevant for studies of 
CP asymmetries in rare hadronic B decays. 

5.4 IR renormalons and power corrections 

The pattern of IR renormalons in the massless limit is very similar to that for finite 
charm-quark mass. Following our analysis in Section ^, we first construct the expansion 
of the distribution function (E^) for small values of r. The result can be written as 



r / 1 

Wl{t, x) = In r + - + 2m 

I — X \ 2 

+ -^ + -(2\ny^-l) + -Aix,T)-{x^{l-x)} 
X X \ X 2) X 
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where 



Mx,t) 



1 ^ 1 , 1 



+ 6{x){ -lnV--lnr + - + - 



(58) 



After integration with the hght-cone distribution amphtude, we find that for the case of 
a symmetric amphtude the leading power-suppressed contributions are given by 
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(59) 



which may be compared with the corresponding resuh in (^). If the distribution amph- 
tude is not symmetric, the function Wl{t) receives a first-order power correction given 
by 



WLiTL) 
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(60) 



As in Section ^, these resuhs could be used to construct a simple model of power 
corrections to factorization in rare hadronic B decays. However, as mentioned earlier, 
this model would be incomplete without taking into account other decay topologies such 
as penguin contractions and interactions with the spectator quark in the B meson. We 
plan to analyze these contributions in a future publication. 



6 Conclusions 

We have used the renormalon calculus to study the asymptotic behavior of the hard- 
scattering kernels entering the QCD factorization formula for the nonleptonic weak de- 
cays B^ D^*'>^M~ . We have obtained explicit results for the Borel transforms and 
momentum distribution functions of the kernels in the approximation of retaining a sin- 
gle renormalon chain (the large- /3o limit). This method estimates power corrections to 
the factorization formula, and allows us to investigate the (divergent) higher-order be- 
havior of the perturbation series for the kernels. From the pattern of singularities in the 
Borel plane, we have derived a simple model of power corrections that is consistent with 
the renormalon analysis. This model accounts only for corrections due to soft, "nonfac- 
torizable" gluon exchange. Other types of power-suppressed effects exist and have been 
estimated in the literature [|I|. 

An unexpected result of our work is that the leading IR renormalon singularity, 
corresponding to a first-order power correction in Aqco/'^^b, vanishes for mesons with a 
symmetric light-cone distribution amplitude. We have explicitly calculated the second- 
order correction and shown it to be numerically small. Higher-order perturbative effects 
are thus expected to be more important than power corrections. We have presented 
analytic results for the contributions of order jSool to the hard-scattering kernels, which 
presumably are the dominant part of the full two-loop contributions. We have also given 
numerical results for the terms of order /^q in the perturbation series. We have shown 
that lowering the renormalization scale below nib improves the rate of approach of the 
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series to their asymptotic values. The BLM scales associated with the imaginary parts of 
the hard-scattering kernels are below 1 GeV (in the MS scheme). This indicates that the 
strong-interaction phases of the decay amplitudes are not insensitive to nonperturbative 
physics. Finally, by showing that the kernels are free of IR divergences and power- 
divergent endpoint singularities, we have proven the factorization formula presented in 
m] to all orders in perturbation theory in the large-/3o limit. 

In the future, it would be worthwhile to carry out a renormalon analysis for the 
phenomenologically more interesting charmless hadronic B decays. Power corrections in 
these processes are, in general, not suppressed by small Wilson coefficient functions. By 
giving explicit results valid in the limit rric — 0, we have accomplished the technically 
most challenging part of such an analysis. It remains to add the contributions from 
penguin contractions and gluon exchange with the spectator quark. 
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Note added: While this paper was in writing, the work Renormalon analysis of heavy- 
light exclusive B decays appeared [|^, in which the authors present expressions for the 



Borel transforms of the hard-scattering kernels that are equivalent to our eq. (|19D , and 
observe the vanishing of first-order power corrections for the case of a symmetric light- 
cone distribution amplitude, in accordance with our eq. (^61) . 
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